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Abstract 

The Lagrangian average (LA) of the ideal fluid equations preserves 
their transport structure. This transport structure is responsible for 
the Kelvin circulation theorem of the LA flow and, hence, for its con- 
vection of potential vorticity and its conservation of helicity. 

Lagrangian averaging also preserves the Euler-Poincare (EP) vari- 
ational framework that implies the LA fluid equations. This is ex- 
pressed in the Lagrangian-averaged Euler-Poincare (LAEP) theorem 
proven here and illustrated for the Lagrangian average Euler (LAE) 
equations. 

PACS numbers: 02.40.-k, ll.lO.Ef, 45.20.Jj, 47.10.-|-g 

Keywords: Fluid dynamics. Variational principles, Lagrangian average 



1 Introduction 

Decomposition of multiscale problems and scale-up In turbulence, 
in climate modeling and in other multiscale fluids problems, a major 
challenge is "scale-up." This is the challenge of deriving models that 
correctly capture the mean, or large scale flow - including the influence 
on it of the rapid, or small scale dynamics. 
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In classical mechanics this sort of problem has been approached 
by choosing a proper "slow + fast" decomposition and deriving evolu- 
tion equations for the slow mean quantities by using, say, the standard 
method of averages. For nondissipative systems in classical mechanics 
that arise from Hamilton's variational principle, the method of aver- 
ages may extend to the averaged Lagrangian method, under certain 
conditions. 

Eulerian vs Lagrangian means In meteorology and oceanography, 
the averaging approach has a venerable history and many facets. 
Often this averaging is applied in the geosciences in combination 

with additional approximations involving force balances (for example, 
geostrophic and hydrostatic balances). It is also sometimes discussed 
as an initialization procedure that seeks a nearby invariant "slow man- 
ifold." Moreover, in meteorology and oceanography, the averaging may 
be performed in either the Eulerian, or the Lagrangian description. The 
relation between averaged quantities in the Eulerian and Lagrangian 
descriptions is one of the classical problems of fluid dynamics. 

Generalized Lagrangian mean (GLM) The GLM equations of An- 
drews and Mclntyre [1978a] systematize the approach to Lagrangian 
fluid modeling by introducing a slow + fast decomposition of the La- 
grangian particle trajectory in general form. In these equations, the 
Lagrangian mean of a fluid quantity evaluated at the mean particle 
position is related to its Eulerian mean, evaluated at the displaced 
fluctuating position. The GLM equations are expressed directly in the 
Eulerian representation. The Lagrangian mean has the advantage of 
preserving the fundamental transport structure of fluid dynamics. In 
particular, the Lagrangian mean commutes with the scalar advection 
operator and it preserves the Kelvin circulation property of the fluid 
motion equation. 

Compatibility of averaging and reduction of Lagrangians for me- 
chanics on Lie groups In making slow -|- fast decompositions and 

constructing averaged Lagrangians for fluid dynamics, care must gen- 
erally be taken to see that the averaging and reduction procedures do 
not interfere with each other. Reduction in the fluid context refers to 
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symmetry reduction of the action principle by the subgroup of the dif- 

feomorphisms that takes the Lagrangian representation to the Eulerian 
representation of the flow field. The theory for this yields the Euler- 
Poincare (EP) equations, see Holm, Marsden and Ratiu [1998a,b] and 
Marsden and Ratiu [1999]. 

Lagrangian averaged Euler-Poincare (LAEP) equations The com- 
patibility requirement between averaging and reduction is handled 
automatically in the Lagrangian averaging (LA) approach. The La- 
grangian mean of the action principle for fluids does not interfere with 
its reduction to the Eulerian representation, since the averaging pro- 
cess is performed at fixed Lagrangian coordinate. Thus, the process 
of taking the Lagrangian mean is compatible with reduction by the 
particle-relabeling group of symmetries for Eulerian fiuid dynamics. 

In this paper, we perform this reduction of the action principle and 
thereby place the LA fiuid equations such as GLM theory into the EP 
framework. In doing this, we demonstrate the variational reduction 
property of the Lagrangian mean. This is encapsulated in the LAEP 
Theorem proven here: 

Theorem 1.1 (Lagrangian Averaged Euler-Poincare Theorem) 

The Lagrangian averaging process preserves the variational structure of 
the Euler-Poincare framework. 

According to this theorem, the Lagrangian mean's preservation of 
the fundamental transport structure of fluid dynamics also extends 
to preserving the EP variational structure of these equations. This 
preservation of structure may be visualized as a cube in Figure 1. As 
we shall explain, the LAEP theorem produces a cube consisting of 
four equivalence relations on each of its left and right faces, and four 
commuting diagrams (one on each of its four remaining faces). 

Euler-Lagrange-Poincare (ELP) cube The front and back faces of 
the ELP cube live in the Eulerian (spatial) and Lagrangian (material) 
pictures of fiuid dynamics, respectively. The top face contains four 
variational principles at its corners and the bottom face contains their 
corresponding equations of motion. The horizontal edges represent 
Lagrangian averaging and are directed from the left to the right. The 
left face contains the four equivalence relations of the Euler-Poincare 



Figure 1: The Averaged EP theorem produces a cube consisting of four equivalence 
relations on each of its left and right faces, and four commuting diagrams (one on each of 
its four remaining faces). 

Theorem on its comers and the right face contains the corresponding 
averaged equivalence relations. Thus, the left and right faces of the 
ELP cube are equivalence relations, and its front, back, top and bottom 
faces are commuting diagrams. 

The back face of the ELP cube displays the LA preservation of vari- 
ational structure in the Lagrangian fluid picture. Hamilton's principle 
with L yields the Euler-Lagrange equations EL in this picture, and 
Lagrangian averaging A preserves this relation. Namely, Hamilton's 
principle with the averaged Lagrangian L yields the averaged Euler- 
Lagrange equations EL. 

This pair of Hamilton's principles and Euler-Lagrange equations 
has its counterpart in the Eulerian picture of fluid dynamics on the 
front face of the ELP cube - whose variational relations are also exactly 
preserved by the LA process. 

The bottom front edge of the cube represents the GLM equations 
of Andrews and Mclntyre [1978a]. Thus, the GLM equations represent 
a foundational result for the present theory. 

The six faces of the ELP cube represent six interlocking equiva- 
lence relations and commutative diagrams that enable modeling and 
Lagrangian averaging to be performed equivalently either at the level 



D. D. Holm 



Lagrangian averaged Euid dynamics 



5 



of the equations, as in Andrews and Mclntyre [1978a], or at the level 
of Hamilton's principle. At the level of Hamilton's principle, powerful 
theorems from other mean field theories are available. An example is 
Noether's theorem, which relates symmetries of Hamilton's principle to 
conservation laws of the equations of motion. Fluid conservation laws 
include mass, momentum and energy, as well as local conservation of 
potential vorticity. The latter yields the Casimirs of the corresponding 
Lie-Poisson Hamiltonian formulation of ideal fluid dynamics and is due 
to the symmetry of relabeling diffeomorphisms admitted by Hamilton's 
principle for fluid dynamics, see Arnold [1966] and Holm, Marsden, 
Ratiu and Weinstein [1985]. In certain cases, the fluid vorticity wind- 
ing number (called helicity - a topological quantity) is also conserved. 
Lagrangian averaging preserves all of these conservation laws.[] Thus, 
the LA Hamilton's principle yields the LA fluid equations in either 
the Lagrangian, or the Eulerian fluid picture, and one may transform 
interchangeably along the edges of the cube in search of physical and 
mathematical insight. 

Remark 1.2 (Eulerian mean) Of course, the preservation of vari- 
ational structure resulting in the interlocking commuting relationships 
and conservation laws of the LAEP Theorem is not possible with the 
Eulerian mean, because the Eulerian mean does not preserve the trans- 
port structure of fluid mechanics. 

Remark 1.3 (Balanced approaches) The LAEP Theorem puts the 
approach using averaged Hamilton's principles and the method of La- 
grangian averaged equations onto equal footing. This is quite a bonus 
for both approaches to modeling fluids. According to the LAEP Theo- 
rem, the averaged Hamilton's principle produces dynamics that is guar- 
anteed to be verified directly by averaging the original equations, and 
the Lagrangian averaged equations inherit the conservation laws that 
are available from the symmetries of Hamilton's principle. 

Outline of the paper In section ^ we begin by briefly reviewing 
the mathematical setting of the EP theorem from Holm, Marsden and 

^We note that the conserved mean topological quantity resulting after Lagrangian 
averaging is the helicity of the mean fluid vorticity, not the mean of the original helicity. 
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Ratiu [1998a,b]. We state the EP theorem and discuss a few of its imph- 
cations for continuum mechanics in vector notation. We also sketch its 
proof, in preparation for proving the corresponding results for the La- 
grangian Averaged Euler-Poincare (LAEP) theorem presented in sec- 
tion ^. Finally, in section ^, we illustrate the LAEP theorem by apply- 
ing it to incompressible ideal fluids. We also mention recent progress 
toward closure of these equations as models of fluid turbulence. 

2 The Euler-Poincare theorem for fluids with ad- 
vected properties 

2.1 Mathematical setting and statement of the EP theorem 

The assumptions of the Euler-Poincare theorem from Holm, Marsden 
and Ratiu [1998a] are briefly listed below. 

• There is a right representation of Lie group G on the vector 
space V and G acts in the natural way on the right on TG x V*: 
{vg, a)h = {vgh, ah). 

• Assume that the function L : TG x V* ^ His right G-invariant. 

• In particular, if oq € V*, define the Lagrangian L^g : TG R 
by Lagi^Vg) = L{vg,ao). Then Lag is right invariant under the lift 
to TG of the right action of Ga^ on G, where Gag is the isotropy 
group of oq. 

• Right G-invariance of L permits one to define i : g x V* —>■ H hj 

i{vgg-^,aog''^) = L{vg,ao). 

Conversely, this relation defines for any i : g x V* ^ H a right 
G-invariant function L : TG x V* ^ H. 

• For a curve git) G G, let 

uit) ^ g{t)9{t)-^ e TG/G = g 

and define the curve a{t) as the unique solution of the linear 
differential equation with time dependent coefficients 

a{t) = -a{t)u{t) (2.1) 
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where the action of u G g on the initial condition a(0) = aQ E V* 



is denoted by concatenation from the right. The solution of (pTl 
can be written as the advective transport relation, 



a{t) = aog{t) 



-1 



Theorem 2.1 (EP Theorem) The following are equivalent: 
i Hamilton's variational principle 

fi2 



6 L^Mt),m)dt = (2.2) 



holds, for variations Sg{t) of g{t) vanishing at the endpoints. 
ii g{t) satisfies the Euler-Lagrange equations for Lag on G. 
ill The constrained variational principle 

5 / i{u{t),a{t))dt = (2.3) 



tl 



holds on g X V* , using variations of the form 



Oti 

6u = — + adur], 6a=-ar], (2.4) 

where r]{t) E g vanishes at the endpoints. 
iv The Euler-Poincare equations hold on gxV* 

2.2 Discussion of the EP equations in vector notation 

When mass is the only advected quantity, the EP motion equation (p.5|) 
and the advection relation (\2-M) for mass conservation are written as 



Here £u denotes the Lie derivative with respect to velocity u, and the 
operations ad* and o are defined using the L2 pairing {f,g) = ffg d^x. 
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The ad* operation is defined as (minus) the L2 dual of the Lie algebra 
operation, ad, for vector fields, namely 

In vector notation ad^ r] is expressed as 

ad„ 1] = rju — UT] = — [u, rj] = — ad,, u = u ■ Vr] — r] ■ Vu . 

The diamond operation o is defined as (minus) the L2 dual of the Lie 
derivative, namely. 



5a 



Here a and 51 /5a are dual tensors and {5l/5a)oa is a one-form density 
(dual to vector fields under L2 pairing). In vector notation the diamond 
operation o for the example of the density D becomes 



Thus, the EP motion equation (|2.5| ) may be written in Cartesian com- 
ponents as 

d 51 d / 5i A 5i du^ ^ / 5i\ _^ 
dt 5u^ J)x^ \5u'^ J 5u^ dx\ dx^ \5D)^ 

and the advection relation for the mass density D eV* satisfies, 

(— + £u) [d <fx^ = , or ^ = - div(Du) . 
Remarks. 

• In passing from coordinate-free forms to their component expres- 
sions, we shall write tensors in a Cartesian basis. For example, 
we shall include the volume form in the mass density by denoting 
it as Dd^x. 
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The EP motion equation and mass advection equation may also 
be written equivalently using Lie derivative notation as 

/ d \ 5i Si / d \ 

The equivalence here of £u and ad* arises because 6i/6u is a one- 
form density and the equality ad*/^ = £ulJ' holds for any one-form 
density [x. 

In the Lie derivative notation, one proves the Kelvin-Noether 
circulation theorem immediately as a corollary, by 

dtJc(u)D5u Jc(u)^dt D 5u J^^^) D 5D 

for any closed curve c{u) that moves with the fluid. In vector 
notation, this is seen as 

d f 1 5£ 
di 



i ^^-dx^/ V;^-cix = 0. (2.6) 
Jc{u) D 5n 5D 



Helicity conservation. In Lie derivative notation, one may 
rewrite the Kelvin circulation theorem as 

d_ 

where 



+ ^y)"^ + dp = 0, 



V — \ • as.= — — • ax , and p = — - — 
D du oD 



Therefore, 

d 



(^— + £u^ {v A dv) = -d (p dv) , 



where A is the exterior product of differential forms. In vector 
notation, this is the helicity equation, 

d 

— (v • curlv) -I- div (u (v • curl v) +p curl v) — . (2.7) 
ot 

Consequently, for homogeneous boundary conditions one finds 

conservation of helicity A = J v Adv = J" v ■ curl v d^x. The helic- 
ity A is a topological quantity that measures the linkage number 
of the lines of curlv, the fluid vortex lines in this case. 
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2.3 Proof of the EP Theorem 

The equivalence of i and ii holds for any configuration manifold and 
so, in particular, it holds in this case. 

The following string of equalities shows that iii is equivalent to iv. 



Finally we show that i and iii are equivalent. First note that the 
G-invariance of L : TG x V* — > R and the definition of a{t) = aQg{t)^^ 
imply that the integrands in ( p^.2D and ( |2.3D are equal. Moreover, all 
variations 5g{t) G TG of g{t) with fixed endpoints induce and are 
induced by variations Su{t) G g of u{t) of the form 6u = drj/dt + ad„?7 
with ri{t) G g vanishing at the endpoints. The relation between 6g{t) 
and ri{t) is given by 7]{t) = 6g{t)g{t)~^ . This is the proof first given in 
Holm, Marsden and Ratiu [1998a]. 



3 Lagrangian averaged Euler-Poincare theory 

We shall place the GLM (Generalized Lagrangian Mean) theory of 
Andrews and Mclntyre [1978a] into the Euler-Poincare framework dis- 
cussed in the previous section. 

3.1 GLM theory from a geometric viewpoint 

The GLM theory of Andrews and Mclntyre [1978a] begins by assuming 
the Lagrange-to-Euler map factorizes as a product of diffeomorphisms. 



Moreover, the first factor g{t) arises from an averaging process. 








(2.8) 



QED 



g{t) = S (t) ■ ~g{t) . 



g{t) = Eit)-git)=git) 



Figure 2: GLM theory factorizes the Lagrange to Euler map at a given time by first 
mapping the reference configuration to the mean position, then mapping that to the 
current position. 



that satisfies the projection property, so that g{t) = g{t). Thus, a 
fluid parcel labeled by xq has current position, 

x^(xo, t) = S it) ■ {g{t) ■ Xq) = S (x(xo, t),t) (current position) , 

and it has mean position, 

x(xo, t) — g{t) ■ xo (mean position) . 

Remark 3.1 Thus, GLM theory first averages the action of the dif- 
feomorphism group, while holding fixed the material objects on which 
the group acts. Then it restores the original action of the group by 
assuming that g{t) ■ g~^{t) — S (i) is also a diffeomorphism. This is 
illustrated in Figure 2. 

The composition of maps g{t) = S (t) ■ g{t) yields via the chain rule 
the following velocity relation. 



^(t)-xo = S(t)-x + rS-(^(t)-xo) 



(3.1) 



Figure 3: The GLM velocities u(x^, t) and u^(x, t) are tangent to the current and mean 
trajectories, and x, respectively. 

By invertibility, xq = g~^{t) ■ = g~^{t) ■ x. Consequently, one may 
define the fluid parcel velocity at the current position in terms of a 
vector field evaluated at the mean position as, 

u(x«, t)=g- g~\t) -^^ = g.g-\t)-^ = u«(x, t) . 

Hence, by using the velocity relation one finds, 

u«(x,t) = S(t)-x + TS-(^r'W-x) 

- f (x,t) + g-u^(x,t). (3.2) 

Here the Lagrangian mean velocity is defined as 

u^(x, t) = u«(x,t) = gg-^{t) ■ x = mg{ty' ■ x = m " xo • (3.3) 

In the third equality we used the projection property of the averaging 
process and found g = g = g from equation (^.11), so that 

u\^,t) = ~g{t)g{t)-' ■^ = u{^,t). 

Thus, the Lagrangian mean velocity vector satisfies = u, so 
is tangent to the mean motion associated with g{t). Hence, one may 
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write equation (|3.2| ) in terms of the mean material time derivative 

D^/Dt as 

u«(x, t) = (I + ■ V) S(x, t)^^E (x, t) . (3.4) 

Likewise, for any fluid quantity x one may define as the composition 
of functions 

X«(x,t) = x(x«,t) = x(2 (x,t),t). 

Taking the mean material time derivative and using the definition of 
D^/Dt in equation (|3.4| ) yields the advective derivative relation, 



As in equation { \i.'S\ ) for the velocity, the Lagrangian mean x^ of any 
other fluid quantity x is defined as 



X^(x, t) = x«(x, t) = x(x«, t) = x( S (t) ■ X, t) = x( ■ xo, t) . 

Taking the Lagrangian mean of equation ( p.5|) and once again using its 
projection property yields 



-y" so that —y'-f^Y 
^ - Dt^ ~\Dt) ' Dt^~\Dt)' 

where x^ = X^~X^ is the Lagrangian disturbance of x satisfying = 0. 

Remark 3.2 T/ie Lagrangian mean commutes with the material 
derivative. Hence, the advective derivative relation (^.d^ ) decomposes 
additively, as 



D. D. Holm 



Lagrangian averaged Buid dynamics 



14 



3.2 Mean advected quantities and their transformations 

Advective transport by g{t) and g{t) is defined by 

a{x.^ ,t) = ao ■ g~^{t) and d{'x, t) = ■ g^^ (t) , 

where = a(xo,0) = a(xo,0), with a,d & V* and the factorization 
g{t) = H (t) • g{t) implies 

a(x, t) = a ■ H (x, t) . 

Note that the right side of this equation is potentiaUy rapidly varying, 
but the left side is a mean advected quantity. 

Since a and d refer to the same initial conditions, Oq, one finds 

ao ■ g^^{t) = a(x, t) = a ■ S (x, t) = a(x^, t) = J?-'(x, t) ■ a^(x, t) , (3.7) 

where J^(x, t) is the tensor transformation factor of a under the 
change of variables S : x ^ x^. For example, one computes formula 
(|3.?1 ) for an advected density as 

(^D(xo)rf^so) ■ = D^{x,t) det{TE) d^x = D{x,t) d^x . (3.8) 

Thus, for an advected density, D, 

D«det(rS)(x,t) = D(x,t), J^(x,t) = det(TS), —D = -d\v{Du). 
For an advected scalar function, s, 

d 

s«(x,t) = 5(x,t) = s^(x,t) , J^=l, —s = -u-Vs. 

For an advected vector field, B, 

K] 5«^ (x, t) = 5^(x, t) , and K'j = det(TS ) (TS -^)} . 
Thus, in the case of an advected vector field, one has 

K-B«(x,t) = B(x,t), with = K = det(TH ) TS -\ 
and an advection relation (e.g., a frozen-in magnetic field) given by 
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Finally, for an advected symmetric tensor 5* one finds 

(T^ ■ "S*^ ■ )ij = Sij , 

whose advection relation is obtained as in the other cases. 

In each case, the corresponding transformation factor JF appears in a 
variational relation for an advected quantity, expressed via equation 



(|3r^ ) as 

Sa^ = 6 {^-^ ■ a) = J^-^ ■ 6a + {6J^-^) ■ a . (3.9) 
This formula will be instrumental in establishing the main result. 

3.3 Lagrangian Averaged Euler-Poincare Theorem 

Let the assumptions hold as listed previously for the EP Theorem 
11 and assume the GLM factorization g{t) = S (t) ■ g{t) with g{t) = 



E{t)-g{t)- Then, 

Theorem 3.3 (LAEP Theorem) The following are equivalent: 
i The averaged Hamilton's principle holds 



6 [\aMt),9{t))dt = (3.10) 



for variations Sg{t) of g{t) vanishing at the endpoints. 
ii The mean Euler-Lagrange equations for La^ are satisfied on G, 



^.TS -^^•TS = (3.11) 
og dt dg 



iii The averaged constrained variational principle 

-12 . 



5 i{ui{t),a^{t))dt = (3.12) 



'ti 

holds on g X V* , using variational relations of the form 

6u^ = TE ■ + adu fj^ + 6E terms , 
6a^ = T^^ • 5a + 5 H terms, and 6a = — afj, 
where fj{t) = 6gg"^ G g vanishes at the endpoints. 
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iv The Lagrangian averaged Euler-Poincare (LAEP) equations hold 
on ^xV* 



Corollary 3.4 (LA Kelvin-Noether Circulation Theorem) 



at D \ou^ ) Jc(u) D \oa^ 
for any closed curve c{u) that moves with the fluid. 



Proof. Via the equivalence of ad* and Lie derivative for a one-form 
density, the (LAEP) equation implies 



for any closed curve c{u) that moves with the fluid. 



3.4 Proof of the LAEP Theorem 



QED 



The equivalence of i and ii holds for any configuration manifold and 
so, in particular, it holds again in this case. To compute the averaged 
Euler-Lagrange equation ( |3.11 ), we use the following variational re- 



lation obtained from the composition of maps g{t) = H (t) ■ g{t), cf 
the velocity relation (^TTj), 

6g{t) = 6E (t) • ~g{t) + TE (t) ■ 6~g{t) . (3.14) 



Hence, we find 



t2 



= 5 LaMt)^9it))dt 



^^^«o Sg+^-^-5g ]dt 



Sg Sg 



' ^_^.TE-^^-^.TE ) -Sgdt. 

og at og 
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This yields the mean Euler-Lagrange equations Here we have 



dropped 5 S— terms, because they do not figure in the variational prin- 
ciple for Lagrangian mean fiuid dynamics at this level of description. 
The following string of equalities shows that iii is equivalent to iv. 



= 6 



i{u^,a^)dt= (^-^ , 6u^^ + (^-^ , 6a^^ dt 



t2 



j^^.TE,^ + ^d,v)-(^.J^-^,£,~a}dt 



In the second line, we again dropped the 5 S— terms and we substituted 
the following variational relations obtained from equations (p.2|) and 



6u^ = TE ■ 6{ijg-^) +6E terms (3.15) 

= TS • + adti f/j + 5 S terms , (3.16) 

6a^ = JF^^- 5a + (52 terms and 5d = —dfj. (3.17) 

Finally we show that i and iii are equivalent. First note that the 
G-invariance of L : TG x 1/* — >■ R and the definition of a{t) = aQg{t)~^ 
imply that the integrands in (|3.10|) and (|3.12|) are equal, both before 



and after averaging. Moreover, all variations Sg{t) G TG of g{t) with 
fixed endpoints induce and are induced by variations 6u(t) G g of u(t) of 
the form 6u = drj/dt + ad^?] with ri{t) G g vanishing at the endpoints. 
The relation between 5g{t) and r]{t) is given by rj{t) = 6g(t)g(t)~^. 
The corresponding statements are also true for the tilde-variables in 
the variational relations ( p.l4| ) and ( p.l5|) - ( p.l7|) that are used in the 



calculation of the other equivalences. 

QED 



Remark 3.5 (Lagrangian Average Conservation Laws/Balances) 

From the viewpoint of the LAEP theorem, the Kelvin circulation theo- 
rem and its associated conservation of potential vorticity for LA flows 
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both emerge because reduction of Hamilton's principle by its relabeling 
symmetries in passing from the material to the spatial picture of contin- 
uum mechanics is compatible with Lagrangian averaging, which takes 
place at fixed fluid labels. 

LA also preserves the kinematic symmetries of Hamilton's principle, 
so, conservation, or balance, laws for momentum and energy for the LA 
dynamics are also guaranteed by Noether's theorem for the averaged 
variational principle, according to its transformations under space and 
time translations. 



4 Application of the LAEP theorem to incom- 
pressible fluids 

4.1 Euler's equation for an incompressible fluid 



For an incompressible fluid, the EP theorem |2 . Ij yields Euler's equations 

as 

for the reduced Lagrangian 

i = j ]^D\u\^ -p{D -l)d^x. (4.2) 

Here the pressure p is a Lagrange multiplier that imposes incompress- 
ibility. The variational derivatives of this Lagrangian are given by 



5i 



j Du-Su+ (^^lu]"^ - p^D - {D - l)Sp d^x . (4.3) 



The expected Euler equation for incompressible fluids is found upon 
setting Z) = 1 in equation ( [4.1|) as 



d 

—u + u ■ Vm + Vp = . (4.4) 

The auxiliary advection relation for the mass density D is the continuity 
equation 

dD 

— = -div(D«), (4.5) 
which, as usual, ensures incompressibility via the constraint D = 1. 



D. D. Holm 



Lagrangian averaged Buid dynamics 



19 



4.2 The Lagrangian averaged Euler (LAE) equations 

The Lagrangian averaged Euler (LAE) equations are derived from the 
LAEP theorem |3.3| as follows. The corresponding averaged Lagrangian 
in the material description is given by 



Do £xo 



|x^P + p« det 



9x« 
9xo 



1 



(4.6) 



Therefore, the reduced averaged Lagrangian in the spatial picture be- 
comes 



(fx 



|u«|2 + p« ( detrS 



D 



(4.7) 



where we have used equation (|3.8| ) in the change of variables. The 
necessary variations of this Lagrangian are given by (dropping the 6 H— 
terms) 



Si 



+ 5p^ ( det TH -D)+ 6p^ detTl 



(4i 



Here we substituted the pressure decomposition p^ = p^ + p^ with 
= p^ and used the projection property of the Lagrangian average. 
Thus, the pressure constraint implies that the mean advected density 
is related to the mean fluid trajectory by 



D = det TS . 

Consequently, the LAE fluid velocity in general has a divergence, 

div u 7^ , 

as was first noticed in Andrews and Mclntyre [1978a]. The Lagrangian 
disturbance of the pressure p^ imposes the constraint 



5/ det TH = . 

This constraint also arises in the self-consistent theory of wave-mean 
flow interaction dynamics in Gjaja and Holm [1996]. It is irrelevant 
here, though, because we are not considering self-consistent fluctuation 
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dynamics. (The self-consistent theory arises from the 5 S— terms that 
we dropped here.) 

The LAE equation may now be written in LAEP form ( |3.13D in 
components as 



and the advected mean mass density D satisfies the corresponding 
mean continuity equation 

dD 

— = -div(Z)u). (4.11) 
When TS = Id + V^, one finds 



v = u5 + ^^^.V^^ = u^-p. (4.12) 

The term p is called the pseudo momentum in Andrews and Mclntyre 
[1978a]. See, e.g., Holm [2001] for a recent discussion and more details. 

Remark 4.1 (Momentum balance) The EP theory of Holm, Mars- 
den and Ratiu [1998a] implies momentum balance in this case in the 
form, 



dt' ' dx^\ ^ V 2 dx 



(4.13) 



exp 



where subscript exp re fers t o th e explicit sp atial dependence arising 
from the S— terms in |u^P = \D^'E/Dt\'^ obtained from equation 



4.3 Recent progress toward closure 

Of course, the LAE equations ( [4.9|) - ( [4.11D are not yet closed. As in- 
dicated in their momentum balance relation (|4.13|) , they depend on the 
unknown Lagrangian statistical properties appearing as the S— terms 
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in the definitions of v and vr. Until tliese properties are modeled or 
prescribed, tfie LAE equations are incomplete. 

Progress in formulating and analyzing a closed system of fluid equa- 
tions related to the LAE equations has recently been made in the 
EP context. These closed model LAE equations were first obtained 
in Holm, Marsden and Ratiu [1998a,b]. For more discussion of this 
type of equation and its recent developments as a turbulence model, 
see papers by Chen et al [1998, 1999a,b,c], Shkoller [1998], Foias et 
al [1999], [2001] and Marsden, Ratiu and Shkoller [2001] and Marsden 
and Shkoller [2001]. An earlier self-consistent variant of the LAE clo- 
sure was also introduced in Gjaja and Holm [1996]. This was further 
developed in Holm [1999,2001]. 

Remark 4.2 (Transport structure) Although the LAE equations 
are not yet closed, their transport structure may still he discussed he- 
cause they are derived in the LAEP context, which preserves the trans- 



port structure. Thus, as in equations ( \2.(^ and (\2. Tj ) we have 



d 



(h V • (ix = , {LAE Kelvin theorem) , (4.14) 

/c(u) 



and 



d f 

— / V ■ curl V d^x = (LAE Helicity conservation) . (4.15) 
at J 

Of course, the LAEP approach is versatile enough to derive LA equa- 
tions for compressihle fluid motion, as well. This was already shown 
in the GLM theory of Andrews and Mclntyre [1978a]. For hrevity, we 
only remark that the LAEP approach also preserves magnetic helicity 
and cross-helicity conservation when applied to magnetohydrodynamics 
(MHD). 
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